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Abstract 

The quantum numbers of the X{3872) meson are determined to be J^'^ = based 
on angular correlations in — )• X(3872)i^'"'" decays, where X(3872) — )• n^TT^J/ip 

^ and J/tp — )• ^^/i^. The data correspond to 1.0 fb~^ of pp collisions collected by the 

LHCb detector. The only alternative assignment allowed by previous measurements, 
J^'-^ = 2 '", is rejected with a confidence level equivalent to more than eight Gaussian 
standard deviations using the likelihood-ratio test in the full angular phase space. 

^ This result favors exotic explanations of the X(3872) state. 
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It has been almost ten years since the narrow X(3872) state was discovered in B'^ decays 
by the Belle experiment Subsequently, its existence has been confirmed by several 
other experiments [2]-^. Recently, its production has been studied at the LHC (5||6]. 
However, the nature of this state remains unclear. Among the open possibilities are 
conventional charmonium and exotic states such as D*^D^ molecules jTj, tetra-quarks fs] 
or their mixtures j^. Determination of the quantum numbers, total angular momentum J, 
parity P, and charge-conjugation C, is important to shed light on this ambiguity. The 
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C-parity of the state is positive since the X(3872) — )■ 7 J/V" decay has been observed 10 

The CDF experiment analyzed three-dimensional (3D) angular correlations in a rel- 
atively high-background sample of 2292 ±113 inclusively-reconstructed X(3872) — )■ 
vr+TT" J/?/' , J/ip — )■ /i"*"/!" decays, dominated by prompt production in pp collisions, 
The unknown polarization of the X(3872) mesons limited the sensitivity of the mea- 



surement of J^*" 12 . A fit of J^^ hypotheses to the binned 3D distribution of 



the J /if) and vrvr helicity angles (6'j/^, 6t,t,) [13 15 , and the angle between their decay 
planes (A0j/^^^ = — ^tttt), excluded all spin-parity assignments except for 1++ or 
2-+. The Belle collaboration observed 173 ± 16 5 ^ X(3872)K {K = or K^), 
X(3872) — > TT+TT^ J/?/^ , J /if) — )■ decays [16]. The reconstruction of the full decay chain 
resulted in a small background and polarized X(3872) mesons, making their helicity angle 
{6x) and orientation of their decay plane sensitive to J^^ as well. By studying 
one-dimensional distributions in three different angles, they concluded that their data 
were equally well described by the l"*"^ and 2 hypotheses. The BaBar experiment 



observed 34 ± 7 X(3872) — uJ/iIj , u — )■ tc^tt'tc^ events 17 . The observed n~^7r~n^ mass 
distribution favored the 2 ^ hypothesis, which had a confidence level (CL) of 68%, over 
the 1"*"^ hypothesis, but the latter was not ruled out (CL = 7%). 

In this Letter, we report the first analysis of the complete five-dimensional angular 
correlations of the B+ X{3872)K+, X(3872) ^ tt+tt- J/ip , J/ip decay chain 

using a/s = 7 TeV pp collision data corresponding to 1.0 fb"^ collected in 2011 by the 
LHCb experiment. The LHCb detector [I8] is a single-arm forward spectrometer covering 
the pseudorapidity range 2 < r/ < 5, designed for the study of particles containing 6 or c 
quarks. The detector includes a high precision tracking system consisting of a silicon-strip 
vertex detector surrounding the pp interaction region, a large-area silicon-strip detector 
located upstream of a dipole magnet with a bending power of about 4Tm, and three 
stations of silicon-strip detectors and straw drift tubes placed downstream. The combined 
tracking system has momentum resolution Ap/p that varies from 0.4% at 5 GeV to 0.6% 
at 100 GeV, and impact parameter (IP) resolution of 20 [xm for tracks with high transverse 
momentum (pt)]^ Charged hadrons are identified using two ring-imaging Cherenkov 
detectors. Photon, electron and hadron candidates are identified by a calorimeter system 
consisting of scintillating-pad and preshower detectors, an electromagnetic calorimeter and 
a hadronic calorimeter. Muons are identified by a system composed of alternating layers 



of iron and multiwire proportional chambers. The trigger 19 consists of a hardware stage. 



based on information from the calorimeter and muon systems, followed by a software stage 



"'^The inclusion of charge-conjugate states is implied in this Letter. 
^We use mass and momentum units in which c = 1. 
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which apphes a full event reconstruction. 

In the offline analysis J/ip — fj^'^fi~ candidates are selected with the following criteria: 
Pt^h) > 0.9 GeV, PTiJ/'ip) > 1-5 GeV, P^r degree of freedom for the two muons 
to form a common vertex, x^^j^(yU+yU~)/ndf < 9, and a mass consistent with the J/ip 
meson. The separation of the J/ip decay vertex from the nearest primary vertex (PV) 
must be at least three standard deviations. Combinations of K'^tt'tt'^ candidates that are 
consistent with originating from a common vertex with x^^.^(i^+7r^7r+)/ndf < 9, with each 
charged hadron (h) separated from all PVs (xip(^) > 9) ^^^'i having pT{h) > 0.25 GeV, 
are selected. The quantity Xip(^) is defined as the difference between the oi the PV 
reconstructed with and without the considered particle. Kaon and pion candidates are 
required to satisfy \n[C{K) / C{7i)] > and < 5, respectively, where £ is the particle 
identification likelihood [20] . If both same-sign hadrons in this combination meet the 
kaon requirement, only the particle with higher is considered a kaon candidate. We 
combine J/ip candidates with K^ti^ti^ candidates to form 5+ candidates, which must 
satisfy Xvtx('V'^-^^^~^^)/^df < 9, pt:{B'^) > 2 GeV and have decay time greater than 
0.25 ps. The J/ip K~^n~7r~^ mass is calculated using the known J/ip mass and the B vertex 
as constraints. 

Four discriminating variables (xj) are used in a likelihood ratio to improve the back- 
ground suppression: the minimal XipW^ Xvtx('^/V'-^^^^^ )/iidf, Xjp(i?+), and the cosine 
of the largest opening angle between the J/ip and the charged- hadron transverse momenta. 
The latter peaks at positive values for the signal as the meson has a high transverse mo- 
mentum. Background events in which particles are combined from two different B decays 
peak at negative values, whilst those due to random combinations of particles are more 
uniformly distributed. The four ID signal probability density functions (PDFs), Vsig{xi), 
are obtained from a simulated sample of B~^ ip{2S)K^ ^ V^(25') -k^-k^ J/ip decays, 
which are kinematically similar to the signal decays. The data sample of B^ — )■ ip{2S)K~^ 
events is used as a control sample for Vsig{xi) and for systematic studies in the angular 
analysis. The background PDFs, V^kgixi), are obtained from the data in the B~^ mass 
sidebands (4.85-5.10 and 5.45-6.50 GeV). We require -2Ylt=iH'Psigixi)/Vhkg{xi)] < 1.0, 
which preserves about 94% of the X(3872) signal events. 

About 38 000 candidates are selected in a ±2cr mass range around the B^ peak in the 
M(J/^ 7r+7r-if+) distribution, with a signal purity of 89%. The AM = M{7r^n- J/tp) - 
M{J/tp) distribution is shown in Fig. [T} Fits to the ipC^S) and X(3872) signals are shown 
in the insets. A Crystal Ball function [2l] with symmetric tails is used for the signal 
shapes. The background is assumed to be linear. The ^(25") fit is performed in the 
539.2-639.2 MeV range leaving all parameters free to vary. It yields 5642 ± 76 signal 
(230 ±21 background) candidates with a AM resolution of a am = 3.99 ± 0.05 MeV, 
corresponding to a signal purity of 99.2% within a ±2.5crAA/ region. When fitting in the 
723-823 MeV range, the signal tail parameters are fixed to the values obtained in the 
ip{2S) fit, which also describe well the simulated X(3872) signal distribution. The fit 
yields 313 ± 26 B~^ — )■ X(3872)i^'"'" (568 ± 31 background) candidates with a resolution of 
5.5 ± 0.5 MeV, corresponding to a signal purity of 68% within a ±2.5crAA/ region. The 
dominant source of background is from B^ — )■ J/tp Ki{1270)~^ , i^ri(1270)+ — t- K^n^-n-' 
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Figure 1: Distribution of AM for i?+ — > J/ip K~^Tr^Tr~ candidates. The fits of the ip{2S) and 
X(3872) signals are displayed. The solid blue, dashed red, and dotted green lines represent the 
total fit, signal component, and background component, respectively. 



decays as found by studying the i^+Tr+vr^ mass distribution. 

The angular correlations in the B'^ decay carry information about the X(3872) 
quantum numbers. To discriminate between the l^"*" and 2~~^ assignments we use the 
likelihood-ratio test, which in general provides the most powerful test between two hy- 
potheses [22j. The PDF for each J^*-^ hypothesis, Jx, is defined in the 5D angular 
space O = (cos6'x, cos ^TTTT, A0x,7r7r, cos ^j/^ , A0x,j/v>) by the normalized product of the 
expected decay matrix element (A^) squared and of the reconstruction efficiency (e), 
V{n\Jx) = \Min\Jx)\^e{n)/I{Jx), where I{Jx) = J |A^(fi| Jx)P e(Q)t/Q. The efficiency 



is averaged over the tt+tt mass (M(7r7r)) using a simulation 23-27) that assumes the 



X(3872) ^ p{770)J/ij, p(770) ^ tt+tt- decay |6||T6)|2|. The observed M{7i7i) distribution 
is in good agreement with this simulation. The lineshape of the p(770) resonance can 
change slightly depending on the spin hypothesis. The effect on e{Q) is found to be very 
small and is neglected. We follow the approach adopted in Ref. [12] to predict the matrix 
elements. The angular correlations are obtained using the helicity formalism, 



\M{n\Jx)\' = E I E 

AA^=-l,+l Aj/^,A^^=-1,0,+1 



'■>X, CJC, 



-(px) X 



D 



where A are particle helicities and D'l_^ are Wigner functions |13 15 . The helicity 
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couplings, Axj^^^x^^, are expressed in terms of the LS couplings |29|,|30], -Bls, where L is 



the orbital angular momentum between the tttc system and the JJip meson, and 5* is the 
sum of their spins. Since the energy release in the X(3872) — )• p{770)J/ip decay is small, 
the lowest value of L is expected to dominate, especially because the next-to-minimal 
value is not allowed by parity conservation. The lowest value for the l^"*" hypothesis is 
L = 0, which implies S = 1. With only one LS amplitude present, the angular distribution 
is completely determined without free parameters. For the 2 ^ hypothesis the lowest 
value is L = 1, which implies S* = 1 or 2. As both LS combinations are possible, the 2 
hypothesis implies two parameters, which are chosen to be the real and imaginary parts 
of a = Bii/{Bii + B12). Since they are related to strong dynamics, they are difficult to 
predict theoretically and are treated as nuisance parameters. 

We define a test statistic t = -2 ln[£(2-+)/£(l++)], where the C{2-+) likelihood is 
maximized with respect to a. The efficiency e{Q) is not determined on an event-by-event 
basis, since it cancels in the likelihood ratio except for the normalization integrals. A 
large sample of simulated events, with uniform angular distributions, passed through 
a full simulation of the detection and the data selection process, is used to carry out 
the integration, I{Jx) Yl^=i \-^{^i\J x)\^ , where A^mc is the number of reconstructed 
simulated events. The background in the data is subtracted in the log-likelihoods using the 
sPlot technique 31 by assigning to each candidate in the fitted AM range an event weight 
{sWeight), Wi, based on its AM value, -21n£(Jx) = -s^2J2ST"u)i \nV{ni\Jx)- Here, 

is a constant scaling factor, = Yl!i=T'' '^i/ Yl!i=T'' '^i^ ■: which accounts for statistical 
fluctuations in the background subtraction. Positive (negative) values of the test statistic 
for the data, tdata, favor the 1"'"+ (2^^) hypothesis. The analysis procedure has been 
extensively tested on simulated samples for the l^"*" and hypotheses with different 
values of a, generated using the EvtGen package f25!|. 

The value of a that minimizes — 21n£(Jx = 2 '",«) in the data is d = (0.671 ± 
0.046, 0.280 ± 0.046). This is compatible with the value reported by Belle, (0.64, 0.27) [16]. 
The value of the test statistic observed in the data is tdata = +99, thus favoring the l"*""*" 
hypothesis. Furthermore, a is consistent with the value of a obtained from fitting a large 
background-free sample of simulated l"*""*" events, (0.650 ± 0.011, 0.294 ± 0.012). The value 
of tdata is compared with the distribution of t in the simulated experiments to determine a 
p-value for the 2 hypothesis via the fraction of simulated experiments yielding a value 
of t > tdata- We simulate 2 million experiments with the value of a, and the number of 
signal and background events, as observed in the data. The background is assumed to 
be saturated by the — > J/ip Ki{127Qi)^ decay, which provides a good description of 
its angular correlations. None of the values of t from the simulated experiments even 
approach tdata; indicating a p- value smaller than 1/(2 x 10^), which corresponds to a 
rejection of the 2 hypothesis with greater than Scr significance. As shown in Fig. § the 
distribution of t is reasonably well approximated by a Gaussian function. Based on the 
mean and r.m.s. spread of the t distribution for the 2 experiments, this hypothesis is 
rejected with a significance of 8.4cr. The deviations of the t distribution from the Gaussian 
function suggest this is a plausible estimate. Using phase-space B^ — )■ J/ip K~^7i~^n~ decays 
as a model for the background events, we obtain a consistent result. The value of tdata 
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falls into the region where the probability density for the 1^"*" simulated experiments is 
high. Integrating the l"*""*" distribution from — oo to tdata gives CL (1"*"^) = 34%. We also 
compare the binned distribution of single-event log-likelihood-ratios with sWeights applied, 
^[7^(^412 ^ , a) /V{fli\l~^~^)], between the data and the simulations. The shape of this 
distribution in data is consistent with the 1"'""'" simulations and inconsistent with the 2 
simulations, as illustrated in Fig. |3j 




100 200 

t = -2 ln[ L{2'yL{r^) ] 

Figure 2: Distribution of the test statistic t for the simulated experiments with J^*^ = 2 ^ and 
a = a (black circles on the left) and with J^'-' = 1++ (red triangles on the right). A Gaussian fit 
to the 2 ^ distribution is overlaid (blue solid line). The value of the test statistic for the data, 
^data; is shown by the solid vertical line. 



We vary the data selection criteria to probe for possible biases from the background 
subtraction and the efficiency corrections. The nominal selection does not bias the M(7r7r) 
distribution. By requiring Q = M{J/iljTnT) — M{J/ip) — M{7ctt) < 0.1 GeV, we reduce the 
background level by a factor of four, while losing only 21% of the signal. The significance 
of the 2 ^ rejection changes very little, in agreement with the simulations. By tightening 
the requirements on the of K and /i candidates, we decrease the signal efficiency by 
about 50% with similar reduction in the background level. In all cases, the significance of 
the 2 ^ rejection is reduced by a factor consistent with the simulations. 

In the analysis we use simulations to calculate the I{Jx) integrals. In an alternative 
approach to the efficiency estimates, we use the i/j[2S)K~^ events observed in the 

data weighted by the inverse of 1 matrix element squared. We obtain a value of tdata 
that corresponds to 8.2a rejection of the 2 hypothesis. 

As an additional goodness-of-fit test for the 1"^"*" hypothesis, we project the data onto 
five ID and ten 2D binned distributions in all five angles and their combinations. They 
are all consistent with the distributions expected for the 1"'""'" hypothesis. Some of them 
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Figure 3: Distribution of —ln[V{fli\2~~^,a)/V{^i\l~^~^)] for the data (points with error bars) 
compared to the distributions for the simulated experiments with J^^ = 1^"*" (red solid histogram) 
and with J^*-^ = 2 ^, a = a (blue dashed histogram) after the background subtraction using 
s Weights. The simulated distributions are normalized to the number of signal candidates observed 
in the data. Bin contents and uncertainties are divided by bin width because of unequal bin sizes. 



are inconsistent with the distributions expected for the (2 a) hypothesis. The most 
significant inconsistency is observed for the 2D projections onto cos6'x vs. cos^jr^. The 
separation between the 1"'""'" and 2~~^ hypotheses increases when using correlations between 
these two angles, as illustrated in Fig. |4j 

In summary, we unambiguously establish that the values of total angular momentum, 
parity and charge- conjugation eigenvalues of the X(3872) state are l"*"^. This is achieved 
through the first analysis of the full five-dimensional angular correlations between final 
state particles in — )■ X (3872) if"*", X(3872) — )• vr+vr" J/-^, J/ip — )■ /i"*"/!" decays using the 
likelihood-ratio test. The 2 ^ hypothesis is excluded with a significance of more than eight 
Gaussian standard deviations. This result rules out the explanation of the X(3872) meson 
as a conventional ric2{l^D2) state. Among the remaining possibilities are the Xci(2^Pi) 
charmonium, disfavored by the value of the X(3872) mass |32j, and unconventional 
explanations such as a D*^D^ molecule jTj, tetraquark state js] or charmonium- molecule 
mixture [O]. 
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Figure 4: Background-subtracted distribution of cos^x for (top) all candidates and for (bottom) 
candidates with | cos6',r7r| > 0.6 for the data (points with error bars) compared to the expected 
distributions for the J^^ = 1^^ (red solid histogram) and J^^ = 2 and a = a hypotheses 
(blue dashed histogram) . The simulated distributions are normalized to the number of signal 
candidates observed in the data across the full phase space. 
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